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ABSTRACT 


A constrained  autoregressive  moving  average  (ARMA)  model 
is  used  to  fit  a train  of  impulses  response  resulting  from 
the  displacement  signal  of  firing  the  M-139  machine  gun.  The 
fitted  constrained  ARMA(2,1)  model  has  only  .38%  residuals 
sum  of  squares  as  compared  with  20.34%  by  the  least  square 
fitting  assuming  a deterministic  model  for  a ten  impulses 
experiment.  The  natural  frequency  of  the  system  can  also 
be  estimated  directly  from  the  fitted  model.  In  a particular 
experiment  where  the  gun  jammed  the  fitted  model  yielded 

a natural  frequency  at  12.6  Hz  while  the  actual  firing  fre- 
quency was  12.58.  Hence  the  fitting  technique  can  be 
employed  to  facilitate  the  design  of  mounting  system. 
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Constrained  ARMA 
Train  of  impulses  response 
Stochastic  differential  equation 
Natural  frequency 
Curve  fitting 
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1 . INTRODUCTION 


To  model  responses  of  a second  order  system  excited 
by  a series  of  impulses,  such  as  the  displacement  of  a 
mounting  system  for  firing  a machine  gun,  one  can  super- 
impose a series  of  impulse  responses  to  fit  a deterministic 
curve.  Assuming  the  errors  being  independent  normal 
random  variables,  the  method  of  least  squares  can  be 
used  to  estimate  the  parameters . However,  in  practice, 
the  errors  are  generally  not  independent;  consequently 
the  estimated  parameters  will  be  inefficient  leading 
to  erroneous  estimation  of  the  parameters. 

The  objective  of  this  paper  is  to  introduce  a time  t \ 

series  technique  which  will  facilitate  the  fitting  of  the 
response  of  a train  of  impulses  data.  The  constrained  autore- 
gressive moving  average  (ARMA)  model  will  be  used  to  fit 
a series  of  impulse  responses,  and  the  mathematical 
justification  will  be  derived.  Since  the  model  contains 
only  the  system  parameters,  the  constrained  ARMA  fitting 

* 4 

is  easier  to  implement  than  the  familiar  deterministic  curve 
fitting  method.  Furthermore,  the  constrained  ARMA  model  can 
be  physically  interpreted.  Real  data  of  the  M-139  machine 
gun  firing  will  be  used  ' o demonstrate  the  constrained  ARMA 
model  fittine  and  its  application  to  solve  a mounting  system 
failure  problem. 
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2.  MODELING  TECHNIQUES 

The  deterministic  model  is  applicable  when  the  curve 
to  be  fitted  is  a known  time  function  containing  a set  of 
unknown  time-invariant  parameters.  The  sampled  values  of 
this  function  at  uniform  time  intervals,  say  A seconds,  are 
assumed  to  be  the  sums  of  the  "true"  values  and  stationary 
white  noise. 

Let  Yfc  denotes  the  N observations  at  discrete  time  index 
t,  and  denotes  a sequence  of  independent  normal  random  vari- 
ables with  mean  zero  and  variance  a 2 (a  ’v/NID  (0  ,a  2),  then 

d t 3 

the  deterministic  model  is  represented  by: 

Yfc  = f (6  , t)  + efc  (1) 

f (0,t)  = the  value  of  the  curve  at  time  t = tA 
- c 

t = continuous  time;  t > 0 

v*  C 

0 = p-vector  of  time-invariant  parameters, 

p > 0 

Note  that  if  f(P,t.  ) denotes  the  continuous  time  function, 
then  f ( 0 , t ) = f(0,tA). 

2 . 1 Response  of  a Second  Order  System  Due  to  a Train 
of  Impulses 

One  frequently  finds  in  practice  that  the  curve  to  be 
modeled  is  a response  of  a second  order  system  excited  by 
an  impulse  or  series  of  impulses.  Let 


3 

C and  u>n  be  the  damping  ratio  and  natural  frequency  of  a 
second  order  system  whose  input  is  a train  of  impulses 

N 

1 cjc<5  (tc_tk)  • where  M denotes  the  total  number  of  impulses 
k=l 

and  ck  and  tk  denote  the  strength  and  time  of  occurrence  of 
the  kth  impulse,  then  the  response  7(8, t)  satisfies  the 
differential  equation 


- — f (8 , t ) + 

2 ~ c 1 


£ ck<S  (tc_tk)  ' 
k=l 

M 

and  f(8,tc)  = l ckG(tc-tk) 

k=l 


(2) 


where  G(tc~tk)  is  the  Green's  function  given  by 


-Co)  (t  -t.  ) sin  u>  /1-c2  (t  -t.)u(t  -t,  ) 
G(Vtk)  =e  2 — - . (3) 


u)n/l-C: 


Here  u(tc~tk)  denotes  the  unit  step  function  whose  value  is 
unity  for  t * tfc  and  zero  otherwise.  An  illustration  of  a 
train  of  impulses  response  is  shown  in  Fig.  1 for  3 impulses. 
Figure  la  shows  the  impulses  of  strengths  c^,  and  c^ 
disturbing  the  system  at  times  t^,  t£  and  t^  respectively. 
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Computed  by  Eq.  9,  Figs,  lb,  lc  and  Id  show  the  respon.  os  of 
the  system  due  to  the  impulses  c^6  (tc~t1)  , c2<S(tc-t2)  and 

(tc~t;j)  in  the  same  order.  Assuming  the  system  is  linear, 
then  the  output  f(B,t)  is  simply  the  sum  of  the  individual 
responses,  as  shown  in  Fig.  le.  Note  that  neither  the  im- 
pulses' strengths  nor  the  time  periods  between  the  impulses 
are  assumed  to  be  uniform. 


2 .2  Deterministic  Curve  Fitting  of  Discrete  Data 

Modeling  of  the  continuous  response  f(B,tJ  given  by 
Eq.  (2)  can  be  conveniently  implemented  if  it  is  observed 
at  a uniform  time  interval.  If  f(S,t) 

denotes  the  observed  values  at  discrete  time  t=l,2,...,N, 
then 

M 

f(&  ,t)  = f(B  ,t  /A)  = Z c.G(tA-t.  ) (4) 

- c Jc=l  k K 

and  the  model  of  Eq.  (1)  can  be  written  as 
M 

Yfc  = l ckG(tA-tk)  + f.fc  (5) 

k 1 

2 

Assuming  that  Et's  are  NIDfO,^  ),  one  can  use  the 
nonlinear  least  squares  method  to  estimate  the  parameters 
C,  wn,  c^  and  t^  for  k=l,2,...,  M. 

As  illustrative  examples,  three  sets  of  response  data 
were  fitted.  It  will  be  discussed  in  Sec.  3 that  the 
digitized  experimental  data  from  firing  the  M-139  machine 
gun  displayed  in  Figs.  2a,  3a  and  4a  are  the  responses  of 
second-order  systems  due  to  1,  3 and  10  impulses  respectively. 


RESIDUALS , SSQ  ■ 20  3<  % 


A 


(d) 

CONSTRAINED  ARMA  FIT  TIN 


6 1 SSQ  = 99  62  % 
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Using  a nonlinear  least  squares  minimization  computer  pro- 
gram to  estimate  the  parameters,  the  best  fit  curves  are 
plotted  in  Figs.  2b,  3b,  4b,  and  three  respective  residuals 
in  Figs.  2c,  3c,  and  4c.  Figures  5a,  5b  and  5c  show  the 
plots  of  estimated  autocorrelation  functions  of  the  residuals 
displayed  in  Figs.  2c,  3c,  and  4c. 


2 . 3 Constrained  ARMA  Fitting  Technique 

It  is  seen  from  Figs.  2a,  2b  and  2c  that  the  deterministic 
model  fits  the  single  impulse  response  data  fairly  well  with 
the  relative  sum  of  squares  of  residuals  (sum  of  squares  of 
residuals  : sum  of  squares  of  the  data  x 100)  of  approximately 
2%.  However,  the  deterministic  model  does  not  satisfactorily 
fit  the  train  of  impulses  response  data  shown  in  Figs.  3a  and 
4a  since  their  respective  residuals  plotted  in  Figs.  3c  and 
4c  have  large  relative  sums  of  squares  (11%  and  20%  respect- 
ively) . Furthermore,  Figs.  5a,  5b  and  5c  shows  that  each 
series  of  residuals  are  highly  correlated,  and  it  appears 
that  the  deterministic  model  is  not  adequate. 

The  dependence  among  et's  can  be  adequately  described 
by  an  ARMA(n,m)  model,  where  n and  m are,  respectively,  the 
autoregressive  and  moving  average  orders  [3].  Hence,  let 
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where 

<}k  ' s = autoregressive  parameters 
9 j ' s = moving  average  parameters 
at's  = orthogonal  decompositions  of  et's. 


a 's  are  NID(0,cf  2) 

t a 

Combining  Eqs.  (5)  and  (6),  it  follows  that 


M 


Y - l c G(tA-t.  ) 

c k=l  K K 


n 


i = l 


M 

l C G(t-lA 
k=l  K 


m 


_tk)1+a 


t“^iejat-j 


(7) 


Since  at's  are  NID(0,oa2),  the  model  given  by  Eq.  (7)  can 
be  fitted  using  the  non  linear  least  square  estimators  in 
the  same  manner  as  the  deterministic  model  given  by  Eq.  (1) . 

It  is  frequently  found  in  practice  that  the  order  of  the 
residuals  et's  can  be  assumed  to  be  two  or  less.  In  fact, 
the  estimated  autocorrelation  functions  plotted  in  Figures 
5a,  5b  and  5c  suggest  that  e^'s  are  of  second -order . There- 
fore, model  (7)  may  be  simplified  by  assuming  that  n=2,  m=l 
and  Eq.  (7)  becomes: 


VjjV V ■ ♦i|Yt-i-kI1ck0<t-u-tk)l 

+ *2(Yt-2")J10kG<t"2''~tk)  1 + at'9lat-l 
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or 


Yt  " *lYf  1 + Vt-2  + at  “Vt-l 


M 

+ l cklG(tA-tk)-*1G(t-lA-tk)-<J>2G(t-2A-tk)]  (8) 

k~  1 


Note  that  fitting  model  (8)  to  the  response  data  is 
rather  inconvenient,  as  there  are  2M+3  parameters  to  be 
estimated.  In  order  to  simplify  this  problem  it  is  assumed 
that  and  G(tA-tk)  are  from  the  same  system,  and  model  (8) 
is  simplified  to  (see  appendix  for  proof) 


Yt  = *lYk-l  + 


^2Yt-2  + at-6lat-l 


(9) 


where  am3  01  are  functi°nally  related  to  z,  and  u>  . 

Let  a = Cwn,  and  b = o^/l-C2,  then 


= 2e  a^  cos  bA 


Q1  = -P+/P2-l, | ex | <1, 

2p  _ b sinh  2aA  -a  sin  2bA 

a sin  bA  cosh  aA-b  sinh  aA  cos  bA 


(10) 

(11) 


(12) 


Although  model  (9)  has  an  appearance  of  the  ARMA  model, 
it  differs  from  the  latter  in  many  respects.  The  most  pro- 
nounced one  is  that  the  parameters  4>2  and  0^  are  func- 

tionally related,  as  clearly  shown  by  Eqs.  (10),  (11)  and 
(12),  while  those  of  the  ARMA  model  are  not.  Furthermore, 


-4 
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Eq.  (11)  shows  that  ip^  in  model  (9)  is  always  negative,  and 
the  value  of  in  the  ARMA  model  does  not  have  this  restric- 
tion. To  distinguish  model  (9)  from  the  general  ARMA(2,1)  model, 
it  is  called  the  constrained  ARMA (2,1)  model  [4]  or  uniformly 
sampled  autoregressive  process  of  order  2 (USA(2))  [3]. 

Note  that  Eqs.  (10,11,  and  12)  show  that  4*^,  $ 2 and  0^ 
can  be  uniquely  determined  from  t, , uj  and  A.  Hence  the  para- 
meter C and  w can  be  estimated  directly  from  digitized  data 
using  model  (9) . 

2 . 4 Examples  of  Constrained  ARMA  Modeling 

The  data  plotted  in  Figs.  2a,  3a  and  4a  was  modeled 
using  Eq.  (9)..  The  least  sciuares  estimators  of  t,  and  u>  are 
obtained  by  Eos.  (9-12),  and  the  continuous 

and  discrete  parameters  are  tabulated  in  Table  1.  The  fitted 
curves  are  plotted  in  Figs.  2d,  3d  and  4d;  the  respective 
at’s  are  plotted  in  Figs.  2e,  3e  and  4e.  Finally,  the  esti- 
mated autocorrelations  of  the  at's  are  shown  in  Figs.  5d, 

5e  and  5f.  Since  the  constrained  ARMA  model  contains  only 
two  parameters,  t,  and  <d  , the  estimation  procedure  is  simple 
and  the  computer  time  required  is  greatly  reduced.  For  ex- 
ample, the  approximated  computer  time  required  to  fit  the 
data  shown  in  Fig.  4a  was  reduced  by  a factor  of  10.  It  is 
also  apparent  from  Figs.  2e,  3e,  4e,  5d , 5e  and  5f  that  the 
constrained  ARMA  fitting  technique  yields  the  curves  that 
fit  the  data  remarkably  well. 

The  constrained  ARMA (2,1)  model  is  used  because  its 
parameters  are  directly  related 


16 

to  the  system  under  studies,  and  these  parameters  can  be 
physically  interpretted . But  if  the  system  parameters  are 
of  secondary  importance,  the  (unconstrained)  ARMA(2,1) 
model  can  also  be  used. 

3.  APPLICATION 

The  three  sets  of  data  used  as  examples  previously  are 
from  experiments  of  firinq  the  M-139  machine  gun  mounted  on 
a table  composed  of  a table  top  and  two  cantilever  beams 
sketched  in  Fig.  fc.  The  mounting  shown  is  adjus- 
table so  that  the  effective  lengths  of  the  beams  can  be 
varied.  The  M-139  machine  gun  was  sot  to  fire  one,  five  and 
ten  rounds,  and  the  table  displacements  were  recorded  during 
these  firing  experiments.  The  data  shown  in  Figs.  2a  and  4a 
are  the  displacement  of  the  table  during  the  experiments  of 
firing  the  M-139  machine  gun  one  and  ten  rounds  respectively. 
There  are  occasions  when  the  M-139  machine  gun  was  set  to 
fire  5 or  10  rounds,  but  it  failed  consistently  after  3 suc- 
cessive fires;  and  the  data  displayed  in  Fig.  3a  is  an  example 
of  the  table  displacements  obtained  from  such  experiments. 

The  displacements  of  the  table  sketched  in  Fig.  6 are 
characterized  by  the  second  order  differential  equation: 


— - ?(8, t )+2tta>  -r—  ?(0,t  )+w  af(B,t  ) = I(t  ) 

dt  2 - c n dtc  ~ c n - c c 


where  I (t  ) is  the  input  force,  and  J, , co  and  T(0,t  ) are 

w II  ^ V, 
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as  described  in  Section  2.  During  the  experiments  of  firing 
the  M-139  machine  gun,  the  input  force  can  be  closely  ap- 
proximated by  a series  of  impulses,  each  with  occurring 
time  t^  and  strength  c^»  k = 1,2,...,M.  It  is  assumed  that 
during  the  course  of  experiments  the  table  vibrates  random- 
ly, hence  the  white  noise,  Z(tc),  is  specified  as  an  addi- 
tional input.  It  has  been  shown  in  Section  2.3  that  the 
discrete  observations  of  the  table  displacement  Yt's  can  be 
modeled  using  Eq.  ( 9 ) , and  the  results  of  constrained  ARMA 
modeling  were  given  in  Section  2.4. 

During  the  course  of  firing  experiments,  the  average 
firing  frequency  of  the  M-139  machine  gun,  defined  as  the 
inverse  of  the  average  time  intervals  between  two  successive 
fires  can  be  estimated.  Figure  7 shows  a plot  of  relative 
displacement  of  the  receiver  and  the  cradle  of  the  M-139 
gun  used  to  estimate  the  average  firing  period,  i.e. 


(2  + inches 


2 secs 
6.5  inches 


.0791 


sec 

fire 


similar  calculations  (using  computer)  were  carried  out  using 
data  from  10  additional  experiments  to  obtain  firing  periods 
of  .0797,  .0809,  .0^97,  .0782,  .0792,  .0809,  .0801,  .0794  and 
.0783  second.  From  these  estimated  firing  periods,  the 
average  firing  frequency  is  12.58  Hz.  For  the  table  mount- 
ings associated  with  the  experiments  in  which  the  M-139  ma- 
chine gun  failed,  the  constrained  ARMA  modeling  technique 


19 

yielded  an  estimated  table's  natural  frequency  of  12.6  Hz 
as  shown  in  Table  1.  Noting  this,  one  can  conclude  that 
during  the  firing  experiments,  the  M-139  machine  gun  failed 
because  its  firing  frequency  and  the  table's  natural  fre- 
quency were  of  the  same  order  of  magnitude.  It  should  also 
be  noted  that  the  table  used  in  the  firing  experiments  which 
the  M-139  machine  gun  did  not  fail  had  different  mountings 
(i.e.,  different  spring  constants).  The  estimated  natural 
frequencies  of  the  table  were  different  from  the  firing  fre- 
quency of  12.58  Hz;  e.g.,  12.0  and  25.7  Hz  as  shown  in 
Table  1. 

5.  CONCLUSION 

1.  It  is  shown  that  a constrained  ARMA(2,1)  model  can 
be  used  to  fit  the  experimental  responses  of  a second- order 
system  due  to  a train  of  impulses  assuming  the  noise  is  also 
an  output  of  the  system. 

2.  The  constrained  ARMA(2,1)  model  is  superior  to  the 
signal  plus  white  noise  model  because  of  (i)  less  parameters, 
(ii)  less  computer  time  and  (iii)  less  residuals. 

3.  The  constrained  ARMA(2,1)  model  can  be  physically 
interpreted  in  terms  of  damping  factor,  c;,  and  natural  fre- 
quency, u>n;  hence  it  can  be  used  in  designing  the  machine  gun 
mounting  system. 
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APPENDIX 

Suppose  that  efc  and  G(tA-tk)  are  from  the  same  system; 
i.e.,  is  a uniform  sampling  of  e(tc)  governed  by  a sto- 
chastic differential  equation 

— — e(t  )+2cu)  -3^-  e (t  ) +u)  2 e (t  ) = Z(t  ) A.l 

j.  2 c n dt  c n c c 

dt  c 

c 

where  Z (tc)  is  the  white  noise  process  with  properties 

E(e(t  ))  = 0,  E (Z ( t ) Z (t  +t ) ) = o 26(t).  The  discrete  pro- 
c c c z 

cess  efc  is  represented  by  [1,2,4], 

Et  = ^let-l  + ^2et-2  + at  ~®lat-l  A. 2 

where  <J>^,<J>2  an<*  ei  are  as  by  Eqs.  (10),  (11)  and  (12). 

Also  G(tA-t^)  defined  by  Eq.  (3)  can  be  written  as 

-a  (tA-t.) 

G(tA-tk)  = e * sin(b(tA-tk))U(tA-tk)/b  A.  3 

Suppose  that  (t-2)A-tk  * 0,  then  U(tA-tk)  = 1,  and 

" 3 A ® 2 A 

using  <J>i  = 2e  cos  bA,  $2  = -e  and  Eq.  A. 3, 
it  can  be  shown  by  algebra  that 

G(tA-tR)  - 4>i  G(t^lA-tk)  - <j>2  G(t^2A-tk)  =0  A. 4 


. . .... 
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Since  G(tA-t^)  = 0 for  tA-t^<0,  it  follows  that 
G (tA-t^) -<t>^G  (t-lA-t^)  — 2^  (t-2A-t^)  = 0 when 

< 0 or  (t-2)A-t^  0.  Hence  Eq.  (8)  can  be  closely 

approximated  by  Eq.  (9);  i.e.. 


Yt  = ^lYt-l  + ^2Yt-2  + at-0lat-’ 


i 
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ABSTRACT 


Force  data  recorded  during  an  experiment  of  firing 
the  M-139  machine  gun  are  analyzed  using  stochastic  diff- 
erential equation  models.  These  models  represent  overall 
M-139  systems,  which  are  composed  of  three  systems:  a 
supporting  table,  a recoil  system  and  a bolt  assembly. 

The  analysis  is  made  to  decompose  these  force  models  to 
obtain  the  parameters  of  the  subsystems. 
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1.  Introduction 

In  a previous  paper  [4],  the  table  displacement  data 
resulting  from  the  experiments  of  firing  the  M-139  machine 
gun  has  been  analyzed  using  a constrained  ARMA(2,1)  model. 
The  model  was  fitted  under  the  assumptions  that  the  sig- 
nals are  a train  of  impulse  responses  and  the  noise  is  due 
to  random  shocks.  The  constrained  ARMA(2,1)  model  fits  the 
table  displacement  data  well,  as  the  residuals  variances 
are  less  than  0.4%  of  the  variances  of  the  data.  The 
damping  ratio  and  natural  frequency  of  the  table  were  es- 
timated from  which  it  was  shown  that  the  M-139  machine  gun 
jammed  when  its  firing  frequency  coincided  with  the  table's 
natural  frequency. 

Physically  a firing  experiment  contains  not  only  the 
table  but  also  the  M-139  machine  gun;  yet,  the  analysis  of 
the  table  displacement  data  did  not  reveal  the  gun  dynamics. 
By  examining  a sketch  of  the  firing  experiment  shown  in 
Fig.  1,  it  can  be  seen  that  the  table  displacement  data 
was  recorded  afar  from  the  gun.  Consequently,  one  expects 
little  or  no  effect  of  the  gun  dynamics  on  the  table  dis- 
placements. During  these  firing  experiments,  the  reaction 
forces  at  the  mounting  of  the  gun  to  the  table  were  also 
recorded.  For  an  experiment  of  firing  ten  sucessive  rounds 
of  the  M-139  machine  gun,  the  recorded  forces,  taken  at  the 
four  corners  of  the  rectangular-shape  gun^s  base,  are  shown 
in  Figs.  (2A-2D) . Since  these  forces  were  recorded  much 
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closer  to  the  gun  than  the  table  displacements,  it  is 
expected  that  these  forces  data  contain  the  gun's  dynamics. 

The  objective  of  this  paper  is  to  fit  the  force  data 
into  linear  stochastic  differential  equation  models.  The 
orders  of  these  models  are  expected  to  be  higher  than  two, 
as  the  force  data  contain  the  gun's  as  well  as  the  table's 
dynamics.  Since  the  differential  equation  fitted  to  the 
force  data  describes  a coupled  system,  a method  of  de- 
composition is  developed  in  order  to  interpret  the  model 
in  the  form  of  the  firing  experiments'  subsystems. 


i 


. 
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2.  Modeling  Technique 

The  digitized  force  data  was  fitted  to  obtain  a 

stochastic  differential  equation  model  which  is  formulated 

by  letting  t denote  a continuous  time,  and  Z (t  ) denote 
c c 

the  white  noise  process  with  properties  E(Z(tc))  = 0 and 

2 

E (Z  (t  ) Z (t  + s)  ) = o 6(s),  where  <5  (s)  is  the  Dirac  delta 

C V Z 

function.  The  continuous  autoregressive  moving  average 

process  of  orders  n and  m (AM(n,m)),  denotes  by  X (t  ),  can 

c 

be  defined  in  terms  of  symbolic  derivatives  as 


d1^  (t_  ) . {?  d"'1 

,.n  c . L.an-i  ..  n-i 
dt  i=l  dt 

c c 


Mtc> 


' 2<tcl+  jSl  bi  Zltc> 


(1) 


where  the  coefficients  a . ' s and  b. ' s are  real. 

Let  y 1 - 1 * be  the  character istic  roots;  i.e.,  p^'s 

axe  solutions  of  the  aquation 


n . n-1 

w +an-llJ  + 


+a,p  tap  - 0 


(2) 


The  required  conditions  for  process  (1)  to  be  stationary  are 
m<n  and  the  real  parts  of  all  p^'s  be  negative  13]. 

The  autocovariance  function  of  X(tc>  is 


y (s)  --  E(X  (tc)X  (tc+s)  ) c j exp  (p  j | s | ) 


(3) 


where 
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cj  = a2z  N(yj)N(-Vj)/  2Uj  H (y*  - u£)(-l)n 


N(u.)  = 1 + b.y  . + b_yf  + 
J 3 ^ J 


k=l,  Mj 


. + b y . 
m j 


When  process  (1)  is  sampled  at  a uniform  discrete  time  inter- 
val, say  A seconds,  the  resultant  discrete  process,  Xfc,  can 
be  represented  by  a uniformly  sampled  autoregressive  moving 
average  model  (USAM)  of  the  form  [2,3,4] 


X - E iX  = a - I 0 a 

^ i=l  1 r 1 t -;  = n 3 t-D 


. 2 
where  at  is  a discrete  white  noise  process  with  variance  a , 

and  4>  ^ ' s and  9.'s  which  are  functionally  related  to  a^’s  and 

bj  ' s defined  by  [2,3] 

n n a n n 

n (l-eMi  B)  = n (1  - A.B)  - 1 - E d>  . B1 
i=l  i=l  1 i=l  1 (6) 


and  if  v_.'s  are  the  n-1  invertible  roots  ( | v_j  | <1 ) of  the 
polynomial 


P(B)  - E c (1-X^ ) n (1-X  B)  (1-A.  B"1) 

j=l  3 3 k=l,k*j  K * 


then 
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n-1  n-1 

n (1-v.B)  = 1 - Z e.B3 

3 j=l  3 (8) 

The  parameters  of  the  continuous  process  (1)  can  be 
estimated  from  a set  of  uniformly  sampled  data  using  max- 
imum likelihood  estimators  or  approximated  by  least  squares 
estimators  [2].  Therefore,  for  a given  set  of  uniformly 
sampled  data,  the  sum  of  squares  of  at's  in  model  (5)  is  ex- 
pressed in  terms  of  continuous  parameters.  Using  this 
procedure,  the  sum  of  squares  of  at' s can  be  minimized  to 
obtain  the  direct  estimates  of  the  continuous  parameters. 
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3.  Modeling  of  Force  Data 

There  are  four  sets  of  measured  forces  in  a firing  ex- 
periment. In  order  to  reduce  these  force  data,  it  is  noted 
from  Fig.  1 that  the  effective  force  is  horizontal,  hence 
the  four  forces  are  averaged  to  obtain  a total  force  series 
as  shown  in  Fig.  2E.  By  examining  this  figure,  one  can  see 
that  the  force  data  should  be  fitted  in  two  separate  aspects? 
that  is,  (i)  force  after  the  firinq  period  and  (ii)  force 
during  the  firing  period. 

The  force  data  were  modeled  using  the  method  summarized 
in  Sec.  2.  The  adequate  model  for  the  force  after  firing 
are  found  to  be  AM  (4,3)  and  for  the  force  during  firing  are 
AM(6,5).  The  continuous  and  discrete  parameters  of  the  two 
models  are  tabulated  in  Table  1.  The  force  after  firing  is 
shown  in  Fig.  3A  and  the  fitted  model  and  at 1 s are  shown  in 
Figs.  3B  and  3C.  Note  chat  the  USAM  (4,3  j models  fits  the 
data  well,  as  the  3um  of  squares  of  at's  is  only  1.7%  of  the 
sum  of  squares  of  the  force  data  after  firing.  The  adequacy 
of  the  fitted  USAM(4,3)  model  can  be  established  from  the 


estimated  autocorrelations  of  the  residuals.  Figure  4A  shows 
that  the  estimated  autocorrelation  values  at  various  lags  are 
below  the  2a  limits  (indicated  by  dotted  lines).  Therefore, 
the  assumption  that  the  at‘ s are  uncorrelated  is  not  violated. 

Figs.  3D-3F  show  the  force  during  firing,  the  fitted 
USAM  (6,5)  model,  and  the  at' s.  The  adequacy  of  the  fitted 
USAM (6,5)  model  is  apparent  from  the  estimated  autocorrelations 
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r 


of  a^.  ’ s shown  in  Fig.  4B.  The  sum  of  squares  of  a^/s  being 
approximately  4.2%  of  the  sum  of  squares  of  the  data. 

It  should  be  noted  that  although  the  models  fitted  to 
the  force  data  are  continuous,  their  values  can  be  estimated 
at  discrete  time  instants  only.  For  this  reason,  the  models 
shown  in  Figs.  4b  and  4e  are  referred  to,  respectively,  as 
U SAM  (4,3)  and  USAM(6,5). 
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4.  Decomposition  of  the  Fitted  AM  Models 

In  order  to  associate  the  fitted  AM  (4,3)  and  AM  (6, 5) 
models  with  the  firing  experiment,  they  are  decomposed  into 
subsystems.  Before  the  decomposition  technique  is  introduced, 
it  should  be  pointed  out  that  when  the  gun  is  operating,  the 
bolt  assembly  moves  relative  to  the  gun.  The  gun  ceases  to 
operate  when  the  bolt  assembly  and  the  gun  are  locked  to- 
gether. The  fitted  models  reveal  this  phenomenon,  as  the 
model  fitted  to  the  force  after  firing  is  of  lower  order 
than  the  model  of  the  force  during  firing. 

4.1  Decomposition  of  the  Force  Model  After  Firing 
(a)  Theoretic  Deterministic  Model 

When,  the  M-139  machine  gun  ceases  to  operate,  the  bolt 
assembly  and  the  rest  of  the  gun  can  be  considered  as  one 
mass.  The  experiment  then  can  be  modeled  as  a two  degree- 
of-freedom  system  as  shown  in  Fig.  5.  Here  M^,  and  k.^ 
represent  the  mass,  damping  coefficient  and  spring  constant 
for  the  table,  and  M2 , C2,  and  k2  represent  the  recoil  sys- 
tem. When  this  system  is  excited  by  forces  F^  and  F^,  the 
system  dynamic,  designated  by  and  X2  can  be  written  in 
terms  of  the  Laplace  transform  variable  as 


For  convenience  let 


then 


T11  - 

M^s2 

+ ic1 

+ C2 ) s + k 

1 + 

T12  = 

C2s 

f k2 

T22 

M2s2 

+ C2 

s + 

k2 

(3) 

'T22F1 

(s) 

+ T12F2 

(s) 

X2  (s) 

,T12F1 

(s) 

* T11F2 

(s) 
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(10) 

(11) 

(12) 


<T11T2  2 


The  measured  force  as  a function  of  time  t is 


x(t)  = C2  -d—  IX2  (t)  - X^t)]  + k2(X2(t) 


T12  > 


2.-1 


(13) 


The  Laplace  transform  of  X(t)  is 
X (s)  = T12  [X2  (s)  - Xl  (s)] 
therefore 


Xl(t”  (14) 


(15) 


% 


■■ 

r3 


(T11T22  - T122)  X<s> 

“ T12  I<T12  * T22)Fl(s)  + (T11  - T12>F2<S» 

In  time  dcmain,  Eq.  (16)  is  of  the  form 


(D4  + a3D3  + a2D2  + OjD  + aQ)  X (t  ) = 

{ k1J«2r2  (t  1 +D  (fcjCj+Cjkj)  F2  (t)  +D2  [ «:1C2+M1k2)  F2  (t)  (t)  1 

+D3[C2(M1F2(t)-M2F1(t)]}  / MxM2  (17) 
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where 


a0 

- kjk./MjM, 

* “l2“22 

(18) 

“1 

- klC2/HlM2 

+k2C1/M1M2 

“ 2CjWjU>22+2 ^2<4)1^u,2 

(19) 

a2 

* k1/M1+k2/M2+k2/M1+C. 

LC2y/MlM2  = 4?1W1C2W2 

(20) 

+o)12+w22  (1+M2/M1) 

“3 

= c1/m1+c2/m2+c2/m1  = 

2c1u1+2c2u>2(1+M2/M1)  , 

(21) 

Cj,  u)^  and  $2»  u2  are'  respectively*  the  damping  ratio  and 
natural  frequencies  of  the  table  and  recoil  system. 

(b)  Comparison  Between  the  theoretic  deterministic 
Model  and  the  fitted  AM (4, 3)  model. 

The  AM (4, 3)  model  fitted  to  the  force  after  firing  is 

(D4  + ct3D3  + a2D2  + ctj^D  + aQ)  X(t) 

= (1  + b^  + b2D2  + b3D3)  Z (t ) (22) 

Note  that  the  right  hand  side  of  Eq.  (22)  cannot  be  interp- 
reted in  terms  of  the  parameters  in  the  right  hand  side  of 
Eq.  (17)  due  to  the  unknown  nature  of  (t)  and  F2  (t)  . 
However*  the  left  hand  side  of  Eq.  (17)  and  Eq.  (22)  are 
the  same;  therefore,  Eqs.  (18-21)  also  define  Oq,  o^*  a2 , and 
c*3  of  Eq.  (22)  in  terms  of  the  M-139  firing  experiment  sub- 
system parameters. 

Since  the  exact  solutions  for  M^,  C^,  k^,  M2*  C2*  and 
k2  using  Eqs.  (18-21)  are  not  possible,  an  approximation 
method  is  developed.  It  can  be  seen  from  Eqs.  (18-21)  that 
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if  and  art  measured,  then  the  values  of  c and 


w2  are  lable.  Since  and  are  not  known,  the 


ratio  M2/M1  is  neglected  using  a priori  knowledge  that  the 


mass  of  the  gun  is  much  smaller  than  the  effective  mass  of 
the  table.  Equations  (20)  and  (21)  can  now  be  written, 


2 2 

<*2  = 4?i<iJi^2tJ2  + Ml  + °°2 


(23) 


2C  ^ co^  "t  2, 2^2 


*24) 


Now  let's  consider  the  autoregressive  roots  y^,  y2,  y^  and 
y4  such  that 


(D-px)  (D-y2)  (D-y3)  (D-y4)  = (DVi^+a^+o^D+ag)  (25) 


It  follows  from  Eq.  (25)  that 


= 


^2^4 


ot,  = 


_^ily2u3  + p 1 lJ  2 v‘  4 + yiy3>-l4  + tJ2y3y4^ 


ou  = 


P 1*J  2 + P 3 + )J  4 + P2|j3  + ^2y4  + p3lJ4 


(26) 


a,  = 


- (Ujl  + P 2 + u 3 + 11 4 5 


: 3 


A 


- 


t 


f 


Comparing  Eqs.  (13,  19,  23  and  24)  with  Sqs.  (26),  it  can 
be  seen  that 

U1/y2  = “?!<*>!  ± (Cj^-l)1^2 

u3'^  ' -s2"2  s “2  <^2-1,1/2  (2 

Therefore,  the  approximate  values  of  the  damping  ratios  and 
natural  frequencies  of  the  table  and  the  recoil  system  can 
be  found  from  the  autoregressive  roots  of  the  fitted  model. 
These  values  are  given  in  Table  2. 
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4.2  Decomposit ion  of  Force  Model  During  Firing. 

The  force  during  firing  is  visualized  to  be  from  a three- 
degree  of  freedom  spring,  mass,  and  damper  system  as  shown  in 
Fig.  6.  Here  k^,  and  Mj  represent  the  table,  k2 , C 2 and 
M2  represent  the  recoil  system  and  k^,  and  represent 
the  bolt  assembly.  Using  a similar  approach  to  estimate  the 
model  it  can  be  shown  that  the  model  of  the  force  during  fir- 
ing is  of  order  six.  Knowing  that  the  mass  of  the  bolt  assemb- 
ly is  much  smaller  than  the  mass  of  the  gun,  the  approximate 
damping  ratios  and  natural  frequencies  of  the  table,  recoil 
system  and  bolt  assembly  can  be  obtained  from  the  autoregress- 
ive roots.  The  values  of  these  parameters  are  also  given 
in  Table  2. 

It  should  be  pointed  out  that  in  a previous  paper  a set 
of  table  displacement  data  from  the  same  firing  experiment 
has  been  analyzed.  The  estimates  of  the  table's  damping  ra- 
tio and  natural  frequency  were  .176  and  25.7  Hz  14].  It  can 
be  seen  from  Table  2 that  the  estimates  of  the  table's  nat- 
ural frequency,  using  models  from  force  signals  after  and 
during  firing  are,  26.4  and  25.9  Hz,  respectively.  Although 
the  force  models  are  decomposed  using  an  approximation  method, 
the  estimates  of  the  table's  natural  frequency  are  well  in 
agreement  with  the  previous  estimate  of  25.7  Hz.  Decomposi- 
tions of  force  models  after  and  during  firing  also  yield 
estimates  of  the  table's  damping  ratio.  The  values  of  these 


estimates,  shown  in  Table  2,  are  . 13  and  .004  as  compared 
with  the  previous  estimate  of  .176. 

5.  Conclusions 

(i)  The  force  data  were  fitted  in  two  parts:  during  and 
after  the  firing  period.  The  adequate  model  for  the 
force  during  firing  is  a continuous  autoregressive 
moving  average  model  of  orders  6 & 5 lAM(6,5)],  with 
the  sum  of  squares  of  residuals  approximately  4.2% 
of  the  total  sum  of  squares.  The  model  fitted  to  the 
force  after  firing  is  AM(4,3)  with  the  residuals  sum 
of  squares  estimated  to  be  1.7%  of  the  total  sum  of 
squares. 

( i i )  The  M-139  experiments  are  visualized  to  be  composed 
of  three  subsystems.  They  are  a table,  a recoil 
system  and  a bolt  assembly.  Assuming  that  the  table's 
mass  is  much  larger  than  the  gun's  mass,  and  that  the 
gun's  mass  is  much  larger  than  the  mass  of  the  bolt 
assembly,  an  analysis  is  made  to  decompose  the  force 
data  models.  As  a result  of  this  decomposition,  the 
damping  ratios  and  natural  frequencies  of  the  three 
subsystems  were  obtained. 

(iii)  The  estimate  of  the  table's  natural  frequency  is 
25.9  Hz  when  the  model  for  the  force  during  the 
firing  period  was  used  and  26.4  Hz  when  the  model 
for  the  force  after  the  firing  period  was  used. 

These  values  are  in  agreement  with  the  estimate  of 
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the  table's  natural  frequency  of  25.7  Hz  using  the 
model  fitted  to  the  table  displacement  data  publisned 
earlier.  Using  the  model  for  the  force  data  durinq 
and  after  firing,  the  table's  damping  ratio  is  esti- 
mated as  .004  and  .13,  respectively.  The  table's 
damping  ratio  obtained  from  the  table  displacement 
model  was  .176. 
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